The Singular Effect of Disorder on Electronic Transport 
in Strong Coupling Electron-Phonon Systems 
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We solve the disordered Holstein model in three dimensions considering the phonon variables to be 
classical. After mapping out the phases of the 'clean' strong coupling problem, we focus on the 
effect of disorder at strong electron-phonon (EP) coupling. The presence of even weak disorder 
(«) enormously enhances the resistivity (p) at T = 0, simultaneously suppressing the density of 
states at the Fermi level, (ii) suppresses the temperature dependent increase of p, and (Hi) leads to 
a regime with dp/dT < 0. We locate the origin of these anomalies in the disorder induced tendency 
towards polaron formation, and the associated suppression in effective carrier density and mobility. 
These results, explicitly at 'metallic' density, are of direct relevance to disordered EP materials like 
covalent semiconductors, the manganites, and to anomalous transport in the A-15 compounds. 



The effect of electron-phonon (EP) interactions in 
weak coupling systems is well understood and can be 
captured within Boltzmann transport theory. However, 
there are materials, e.g, the A-15 compounds, where the 
electron-phonon (EP) coupling is large and the interplay 
of disorder and EP interaction leads to several anoma- 
lous features [1-4] in the resistivity, p, including a rapid 
rise in the residual resistivity at weak disorder, a sharp 
reduction in the density of states (DOS) at the Fermi 
level, suppression of the temperature (T) dependence of 
p and, sometimes, a regime with dp/dT < 0. The con- 
junction of EP coupling and disorder has dramatic effects 
on transport in covalent semiconductors as well [5], and 
also controls the resistivity and optical spectra in the low 
T spin polarised phase of the manganites [6] . 

These transport anomalies, ranging across a wide va- 
riety of materials, have a common origin in the interplay 
of polaronic tendency (arising from strong EP coupling) 
with quenched disorder and thermal fluctuations. How- 
ever, to capture these effects within a unified framework 
we need to handle the strong coupling and quenched dis- 
order non perturbatively, and simultaneously retain the 
spatial correlations in the problem. 

We accomplish that in this paper, solving the Holstein 
model in three dimensions, with arbitrary coupling and 
quenched disorder, through a new real space Monte Carlo 
(MC) technique [7] explicitly retaining the spatial corre- 
lations and localisation effects in the problem. We work 
with adiabatic phonons, i.e, treat the phonons as an- 
nealed classical variables. Apart from providing the first 
complete solution to this problem, our results specifically 
highlight: (i) the huge amplification of quenched disor- 
der by strong EP coupling, and the dramatic increase 
in residual resistivity with weak disorder, (ii) the rapid 
suppression of the 'thermal' component of p{T) with in- 
creasing disorder, and (Hi) the emergence of a regime 
with dp/dT < 0. We correlate these effects with (a) the 



appearance of a weak pseudogap in the density of states 
and (b) the suppression of the optical spectral weight, 
suggest a transport phenomenology for these disordered 
EP systems, and compare our results with data on the 
disordered A-15 compounds. 

Model: The disordered Holstein (d-H) model with spin- 
less fermions and classical phonons is described by: 



(ij) 



Z^ 



{ei- p)ni- \^niXi+ Hk (1) 



The t are nearest neighbour hopping on a simple cubic 
lattice, Ei is the quenched binary disorder, with value 
±A, II is the chemical potential, A is the EP interaction, 
coupling electron density Ui to the local distortion Xi, and 
the 'restoring force' arises from Hk = {K/2) J^i^f- The 
parameters in the problem are X/t, A/i, electron density 
n, and temperature T. We measure energy, frequency, 
T, etc, in units of t, and use K — 1. 

While the strong coupling problem with one electron 
has been widely studied, we do not know of controlled 
approximations to handle strong coupling at metallic 
densities [8], except within dynamical mean field theory 
(DMFT). DMFT has been used in the adiabatic hmit 
[9] to map out the Fermi liquid and polaronic insulator 
regime, and the occurence of charge order (CO) [10]. The 
effect of strong EP coupling on p{T), including the effect 
of disorder, has also been explored [11]. 

Our present understanding of strong coupling EP sys- 
tems owe much to DMFT, but this approach has lim- 
itations in the presence of disorder, (i) Disorder and 
strong EP interactions reinforce each other and generate 
a strong local correlation between e^ and Xi, inaccessible 
in the statistically homogeneous description of DMFT. 
(ii) In such a system the conductivity can no longer be 
computed in terms of the single particle self energy, as 
in DMFT, and vertex corrections leading to localisation 
are crucial to recover the dp/dT < regime. 



With the possibiUty of polaron formation, and the need 
to handle disorder and thermal fluctuations, MC calcu- 
lation is the only unbiased method for approaching the 
d-H problem. However, the cost of annealing the phonon 
variables, via the standard exact diagonalisation based 
MC (ED-MC), grows rapidly with system size. For a 
cube with N = L^, the ED-MC cost grows as N'^, allow- 
ing typically L — A with current resources. We instead 
use a "travelling cluster approximation" (TCA) [7] for 
annealing the phonons. Since TCA avoids iterative di- 
agonalisation of the electron Hamiltonian, and estimates 
the energy cost of MC moves based on a smaller 'cluster' 
Hamiltonian, we can access sizes upto 12^ with relative 
ease. We use a 4^ cluster for the TCA based updates. 
Once the phonon variables are equilibriated, we diago- 
nalise the full electron Hamiltonian in the equilibrium 
configurations to compute [12] electronic properties. We 
use A^ — 8'^, checking with 10^ for size effects. Our re- 
sistivity is in units of ?iao/(7re^), where oq is the lattice 
spacing. As a rough measure, p ^ 100 in our units is 
equivalent to the Mott resistivity, about 2 — 3 mr^cni. 

Let us start with the phases. In the absence of dis- 
order, the finite T phases can be broadly classified as 
(i) Fermi liquid (FL) with dp/dT > 0, (ii) a polaronic 
liquid (PL), with 'large' lattice distortions, no long range 
positional order, and dp/dT < 0, and (Hi) the charge or- 
dered insulator (COI) state. Unlike at T = [13], all the 
phases have finite density of states at the Fermi level. 
There is no 'phase boundary' between FL and PL and 
the distinction is only in terms of dp/dT. 

For A < 2.0, Fig.l.(a), the only feature in the n — T 
phase diagram is the occurence of CO at n = 0.5, and 
the coexistence region separating the CO phase from the 
FL. With increasing A the CO region extends down to 
n ~ 0.35 (a chessboard pattern with 'vacancies' [13]), 
while Tco aX n = 0.5 grows and then diminishes slowly 
beyond A ~ 3.5. We will discuss the CO phase in detail 
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FIG. 1. Phase diagram: clean model. Panels (a) — (d) . show 
the n~T phase diagram with increasing A. The "chessboard" 
pattern represents a {7r,7r,7r} CO phase, and the hatched re- 
gions indicate phase separation. 



elsewhere [15]. The coexistence width at T = between 
the FL and COI increases with increasing A but the tem- 
perature window for coexistence reduces since the finite 
T strong coupling FL itself has strong density fluctua- 
tions. By the time A = 3.0 the FL regime has almost 
vanished, the PL and COI phases dominate the phase 
diagram, and the system is insulating at all n and T. 

The interesting regime for studying the impact of disor- 
der and thermal fluctuations is roughly A ^ 2 — 3 where 
the system is at strong coupling but the metallic state 
still survives. Our studies will be focused on A = 2.0, 
strong coupling but still metallic, and n — 0.3, and we 
will use A = 1.0 as our 'weak coupling' reference. To set 
the scale, the single polaron threshold is [13,14] Ac = 3.3, 
while at n = 0.3 collective polaronic localisation [13] sets 
in at A « 2.5. We study the impact of disorder staying 
below this threshold. 

The presence of disorder in the FL phase tends to cre- 
ate density inhomogeneities, Srir- This effect is weak 
at weak disorder in an uncorrelated system but can be 
strongly amplified by a positive feedback if the EP cou- 
pling is large. The phenomenon is understood in the 
context of a single electron [5], but an understanding of 
the dense many electron state, and transport properties, 
has remained out of reach and is our main focus. 

The resistivity arises from a combination of quenched 
disorder and lattice distortions. For a phonon config- 
uration {a;}", the "potential" seen by the electrons is 
^f — ei — \xf. Subtracting out the spatial average, 
fct, the fiuctuating part is (5^f = rjf ~ ei — \{xf — Xa)- 
Since the spatially averaged distortion does not depend 
on configurations, i.e, Xa ~ x, we can write xf — x^ ~ 



(^° 



{xf — a;?), where x^ is the T — distortion at 



Ri. Let us define 6x^ 



X and dxf 



so that (5x° is the T = distortion at a site (with 
respect to the spatial average), and Sxf is the ther- 
mal fiuctuation about that local distortion. In terms of 
these, the fiuctuating background seen by the electrons 
is: rjf r^ ei ~ A((5x^ -I- Sxf). Averaging spatially, and over 
equilibrium configurations, the variance of this effective 
'disorder' is ry^ = ((e, - AfeO)^) + X^{{Sxf)^) - 2A((e, - 
XSx'i)Sxf) = fJo + T]'^ + Vcorri ^ay. ?7g is a measure of 
disorder at T = 0, ryj is a rough measure of "thermal 
disorder" , and rj'^g„ locally correlates the thermal fluctu- 
ations and the T = disorder. 

In the clean FL there are no lattice distortions at 
T = 0, i.e, feo = 0, so 772 = A2((fef)2). The stiflF- 
ness to thermal fluctuations of the Xi, at low T and weak 
coupling, can be roughly estimated from the energy cost 
{K/2)xf + X'^x{''T')x'i, where x is the local, zero frequency, 
density response function of the electron system. The ef- 
fective stiffness is Kgff — K + 2A2x(n) and ({Sxf)'^) oc 
T/Keff. The thermal disorder is r/^ = rj^ ^ X^T/K^ff. 
Electron scattering from these fluctuations would lead to 
p{T) ex r^f,, for rjT ^ W, where W = 12t is the bare 
bandwidth, but as rj^ grows the resistivity rises faster 



than Tj^, due to localisation corrections, unlike in DMFT. 
The nature of ry^ and p{T) in the clean limit and A = 1.0, 
and approximately at A = 2.0, Fig. 2, can be understood 
within the above scenario. 

For A 7^ the effective disorder at T = arises from 
ei and the induced lattice distortion: rj'^ = {{a — XSx'^)'^) . 
If Sx'^ is small, as one expects at weak disorder and weak 
coupling, then r]"^ ss (e|) = A^, leading to p(Q) ~ A^. 
We have checked this at A = 1 and n = 0.3, Fig. 2. (a) in- 
set, and the expected dependence clearly holds. However, 
when A = 2, even weak disorder can create strong density 
inhomogeneities and induce large quenched distortions 
[5]. In this regime, 77^ is dominated by A^(((5x°)^). Since 
Sx° ~ iX/K)Sn°, we have t?^ ^ {X^ / K'^){5r4f . Monte 
Carlo annealing of the phonon variables reveals that den- 
sity inhomogeneities can be strong, 5n^ ^ C(l), and the 
A* factor leads to a huge amplification of the external 
disorder. The rapid growth in /o(0), Fig.2.(c) inset, arises 
due to this. The small deviation from p{Q) ex A^, at 
small A^, is a finite size effect. 

The effect of disorder and thermal fluctuations is ad- 
ditive in both rf and p{T) at A = 1.0, as borne out by 
the parallel curves in Fig. 2. (a) and Fig.2.(c), i.e, p{T) 
obeys Mathiessens rule. However at A = 2.0 the same 
A leads to much siron^er effective disorder, rf'(T), com- 
pared to A = 1.0. What is more striking is that despite 
the drf /dT > in both cases the "thermal increase" in 
p{T) at A = 2.0 is actually suppressed with increasing A, 
Fig.2.((i), leading eventually to a regime with dp/dT < 0, 
by the time A ~ 0.5. We will argue that these transport 
anomalies at A = 2.0 arise from the disorder induced po- 
laron formation tendency, at T = 0, and its weakening 
with increasing temperature. We next discuss 
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FIG. 2. (Colour online) Effective disorder, i]^ , and resis- 
tivity, p(r), at n = 0.3. (a) - (b). shows r;^: (a). A = 1.0, 
(6). A = 2.0. (c) — (d). p{T), corresponding to the respective 
panel above. Inset to (a). ?7^(A) at T = 0, black A = 1.0, red 
A = 2.0, inset to (c). p(0, A), black A = 1.0, red A = 2.0. 
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FIG. 3. (Colour online) Density of states, at A = 2.0, 
n = 0.3. The A and T are marked on the panels, (c). Temper- 
ature dependence of A''(0), the DOS at ep, (d). Correlation, 
at T = 0, between N{0) and p(0), varying A. 

the DOS and optics, and then our overall scenario. 

The DOS, N{u!, T), is featureless in the clean problem. 
Fig. 3. (a), even at A = 2.0, typical of gradual thermal 
disordering of a FL. The presence of weak disorder, e.g, 
A — 0.6, Fig. 3. (6), leads to the formation of a weak pseu- 
dogap around uj = 0. The transfer of spectral weight to 
lower frequencies in N{u!) suggests possible localisation 
of some electrons into polaronic states. The pseudogap 
fills quickly with increasing T and, for A ^ 0.8, vanishes 
for T > 0.15, Fig.3.(c). The suppression of DOS at ep 
and the increase in residual resistivity, with increasing 
disorder, have a monotonic relation, Fig.3.((i). 

The optical conductivity, cr(a;), in the clean problem. 
Fig. 4. (a), is Drude like as expected of a Fermi liquid. 
Very weak disorder, A = 0.3, Fig.4.(&), retains the Drude 
character, but suppresses the overall magnitude. At 
A = 0.6, however, where anomalies were visible in p{T) 
and N{ui), the optical response is non Drude at all T, 
suggesting the existence of localised states. The "effec- 
tive carrier number" neff{uJ,T) — L a(Lu,T)dLo, which 
is a rough measure of the kinetic energy, is reduced by a 
factor ~ 10 (at (I; = 1) as A increases from zero to 0.6 
at T = 0, Fig.4.(d). With increasing T the strongly lo- 
calised particles delocalise, as borne out by the increasing 
ne//(T), Fig.4.(d)., and the DOS at ep, Fig.3.(6). 

We propose the following approximate "two fluid" 
framework to approach the results of Figs. 2-4. (a). At 
T — the clean adiabatic EP problem allows only two 
kinds of collective states [13] at n = 0.3, (i) a band metal 
(with no lattice distortions) for A ^ 2.5, or (ii) a pola- 
ronic insulator, with all electrons localised, for A ^ 2.5. 
The presence of weak disorder at A = 2.0 induces strong 
localisation of a fraction, fp^i = npoi/n of charge carriers. 
npo;(A) « J^ P{x, A)dx, where P{x, A) is the 
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FIG. 4. Optical response; (a) — (c): a{u>,T) at A = 2.0 and 
n = 0.3. (a) A = 0, (6) A = 0.3, (c) A = 0.6. Panel (d) 
T dependence of rieff, the integrated low frequency optical 
spectral weight (see text) for various A. 

distribution of lattice distortions and Xmin is the upper 
edge of the clean FL peak. However, states near ep con- 
tinue to be delocalised. The fraction fpoi increases with 
increasing disorder [15], suppressing rie// and A^(0) and 
increasing p{0). The fraction, 1 — fpoi, of extended states 
are themselves strongly scattered by the large distortions, 
Xi, but still maintain 'metallic' conduction, (b). With 
increasing T, the electrons in extended states scatter 
off thermal fluctuations in Xi and contribute to a grow- 
ing p{T), while localised states close below ep provide 
a "parallel" conduction channel, whose conductivity in- 
creases with increasing T. The reduction in dp/ dT at 
intermediate A comes from these competing tendencies. 
For A > 0.8 there are few extended states and at low 
T the delocalising trend dominates. We have explicitly 
checked the T dependence of rir, [15], and the density 
inhomogeneity, due to strongly localised electrons, weak- 
ens with increasing T . At higher T, after the polaronic 
effects have disappeared, p{T) again roughly tracks rf'. 

Some care is needed in applying our single band adia- 
batic results to real materials, where (?) both the band- 
structure and the EP coupling could be more compli- 
cated, and (m) the phonon frequency, Uph, although usu- 
ally ^ t is still finite. While the former will affect any 
quantitative comparison, the later brings in the phys- 
ical effect of phonon assisted hopping, and possibly a 
larger conductivity from the 'localised' states than in the 
^ph = limit. Such a {iOpu ^ 0) mechanism had been 
studied in the limit of strong disorder and weak coupling 
[16] to explore the sign change in dp/ dT. Our focus is on 
the complementary strong coupling, weak disorder limit. 

The interplay of disorder and strong EP coupling 
is best documented in the 'high Tc' A-15 compounds 
[1-3,17]. There, in addition to 'resistivity saturation' [18] 
in the clean limit, the combination of strong EP coupling 
with weak disorder leads to the following effects: {i) p(0) 
increases enormously^ in response to disorder (alpha par- 
ticle damage) [1] compared to similar disorder in a 'weak 



couphng' material, Nb, say, [ii) as p(0) grows, the ther- 
mal increase in p{T) is systematically reduced, tending 
to a limit dp/dT ^^ at strong disorder [1], [Hi) there is 
a degradation of Tc with increasing /9(0), suggesting a re- 
duction in the DOS at ep [2,3]. {iv) Some other systems, 
e.(7, LuRh4B4 actually go over to dp/dT < with increas- 
ing p{0) [4]. All the features, (i) — (iv) above, find a con- 
sistent explanation in terms of our results. The reduction 
in DOS inferred from Tc degradation bears a close paral- 
lel to our correlation between A^(0) and p(0), Fig.3.(d). 

Conclusion: In this paper we solved the problem 
of transport in a dense, strong coupling, disordered 
electron-phonon system, explicitly in three dimensions. 
Our results reveal how disorder affects the many elec- 
tron system, close to a collective polaronic instability, by 
strongly localising a fraction of electronic states, and dra- 
matically reducing the mobility of extended states that 
survive near the Fermi level. We mapped out the inter- 
play of these polaronic and extended states on transport, 
spectral, and optical properties, suggested a phenomenol- 
ogy for these class of systems, and compared our results 
with data on the A-15 compounds. 
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